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Abstract 

The  charge  characteristics  of  a  lead- acid  cell  are  observed  by  using  a  numerical  simulation.  The  effect  of  various  parameters,  such  as 
the  concentration  exponent  of  charge  reaction,  morphology  parameter  and  limiting  current  density,  on  the  cell  voltage  during  charge  has 
been  investigated  by  including  the  dissolution-precipitation  mechanism  of  the  negative  electrode.  As  the  charging  current  density  is 
increased,  the  concentration  gradient  is  increased  due  to  the  high  resistance  of  the  electrolyte  migration,  especially  at  the  interface 
between  the  positive  electrode  and  the  reservoir.  ©  2000  Elsevier  Science  S.A.  All  rights  reserved. 
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1.  Introduction 

One  of  the  important  factors  needed  to  use  the  batteries 
for  electric  vehicles  is  the  possibility  of  fast  charging.  For 
fast  recharging,  high  currents  are  required.  This  method, 
however,  is  expected  to  cause  excessive  heating  of  the 
batteries.  Overheating  severely  reduces  cycle-life.  There¬ 
fore,  a  method  is  required  to  charge  the  batteries  quickly 
without  damage. 

Maja  et  al.  [1]  examined  the  effect  of  different  parame¬ 
ters,  such  as  plate  dimensions,  amount  of  acid,  porosity  of 
active  material  and  modality  of  charging,  on  the  charge 
performance.  Valeriote  and  Jochim  [2]  studied  the  key 
components  to  reduce  temperature  increase  for  charging 
the  battery.  They  found  that  the  temperature  rise  was 
determined  not  only  by  the  amount  of  heat  produced  by 
the  charge  process  but  also  by  the  rate  of  heat  dissipation 
and  the  heat  capacity  of  the  battery.  The  effects  of  charg¬ 
ing  current,  depth-of-discharge  and  battery  design  were 
observed  in  the  charge  process.  Chang  et  al.  [3]  discussed 
fast-charging  effects  on  hybrid  lead-acid  batteries  which 
contained  antimonial  positive  grids  and  non-antimonial 
negative  grids.  They  used  the  constant  resistance-free  volt¬ 
age  which  compensated  for  the  ohmic  voltage  drop  to 
reduce  the  internal  resistance  of  the  battery. 
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Gu  et  al.  [4]  compared  the  voltage  penalty  and  time  to 
the  voltage  cut-off  at  moderate  and  low  temperature,  re¬ 
spectively.  Gu  et  al.  [5]  developed  the  model  not  only  to 
account  for  coupled  processes  of  electrochemical  kinetics 
and  mass  transport  occurring  in  a  battery  cell,  but  also  to 
consider  free  convection  resulting  from  density  variations 
due  to  acid  stratification.  These  two  investigations  only 
considered  the  charge  transfer  in  the  electrode  kinetics  of 
the  negative  electrode  during  charge  and  did  not  take  into 
account  the  solid-state  reaction. 

Ekdunge  et  al.  [6]  showed  that  the  cathodic  polarization 
curves  at  different  states  of  charge  exhibited  limiting-cur- 
rent  phenomena  due  to  rate-determining  dissolution  of  lead 
sulfate  or/and  diffusion  of  lead  ions. 

In  this  study,  we  have  incorporated  the  solid-state  reac¬ 
tion,  that  is,  the  dissolution  of  lead  sulfate  and  the  diffu¬ 
sion  of  lead  ions,  and  have  developed  a  mathematical 
model  to  predict  accurately  the  effect  of  various  electrode 
parameters  on  the  charge  performance  of  a  lead- acid  cell. 


2.  Mathematical  model 

2.7.  Governing  equations 

The  approach  presented  here  was  based  on  a  macro¬ 
scopic  model  [4,7]  and  a  mathematical  model  [8,9]  which 
describe  the  discharge  performance  of  a  lead-acid  cell. 
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The  assumptions  made  in  the  development  of  the  model 
are  follows. 

(i)  The  lead-acid  cell  consists  of  a  lead  dioxide  elec¬ 
trode  (Pb02),  an  electrolyte  reservoir,  a  porous  separa¬ 
tor,  and  a  lead  electrode  (Pb). 

(ii)  Each  electrode  is  a  porous  and  electronically  con¬ 
ductive  matrix  with  pores  occupied  by  sulfuric  acid 
solution. 

(iii)  The  model  is  one-dimensional  in  the  direction 
perpendicular  to  the  plane  of  the  electrodes. 

(iv)  Porous  electrodes  are  assumed  to  be  macro-homoge¬ 
neous. 

(v)  The  cell  can  be  considered  to  be  isothermal  during 
its  operation. 

During  charge,  the  electrode  reactions  in  a  lead-acid 
cell  are  expressed  as  follows: 

PbS04(s)  +  2H20  ->  Pb02(s)  +  HS04“  +  3H+  +  2e“ 

(i) 

for  the  positive  electrode,  and 

PbS04(s)  +  2e“  +  H+  ->  Pb(s)  +  HS04“  (II) 

for  the  negative  electrode.  Lead  sulfate  is  converted  to  lead 
dioxide  (positive  electrode)  and  lead  (negative  electrode) 
during  charge. 

The  kinetic  behaviour  of  a  lead  electrode  in  a  lead- acid 
cell  is  influenced  by  local  activity,  overpotential,  elec¬ 
trolyte  concentration,  and  current  density  distribution  [6]. 
The  reaction  kinetics  of  the  negative  electrode  assume  that 
the  charge  reaction  consists  of  two  or  more  elementary 
reactions:  dissolution  of  lead  sulfate,  diffusion  of  lead  ions 
to  an  active  lead  surface,  and  electrochemical  reactions. 
The  transfer  current,  which  the  solid-state  reaction  is  in¬ 
cluded  in,  is  written  as  follows  [8,9]: 


The  current  density  in  the  solid  phase,  ,  follows  Ohm’s 
law.  The  current  density  in  the  electrolytic  solution,  i2,  is 
proportional  not  only  to  the  concentration  gradient  but  also 
to  the  electric  potential  gradient. 

The  material  balance  for  the  acid  concentration  in  the 
liquid  phase  is  given  by  convection,  diffusion  and  migra¬ 
tion  of  mobile  ionic  species.  The  volume-average  velocity 
v  *  is  used  as  the  reference  velocity  [10].  The  effective 
diffusivity  D  '  is  given  by  the  diffusion  coefficient  of  the 
electrolyte,  porosity,  and  tortuosity  [11].  The  porosity  of 
the  electrode  during  charge  is  changed  by  the  electrode 
reaction  and  is  proportioned  to  the  transfer  current. 

The  overall  electrode  reaction-rate,  j,  in  the  positive 
electrode  is  represented  by  the  Butler- Volmer  equation. 
The  overpotential  r)  at  this  electrode  is  defined  as: 

^7  =  01  02  kpb02  (3) 


where  <pl  is  the  potential  of  the  solid  phase  and  c/>2  is  the 
potential  of  the  liquid  phase.  UPh0i  denotes  equilibrium 
potential  evaluated  at  a  reference  concentration  cref. 

The  electroactive  surface  area,  a ,  can  be  related  to  the 
state-of-charge.  Note  that  the  charging  reaction  should  stop 
when  lead  sulfate  formed  on  the  previous  step  is  com¬ 
pletely  converted  [4],  i.e., 

/  n  \  l  q  \ ?  /  n  \ 


a 


1  - 


Q 


Q 


=  a 


max 


max  / 


Q 


x 


\  ^  max  / 


l  - 


Q 


Q 


(4) 


max  / 


where  am„Y  denotes  the  maximum  active  surface  area  of 
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the  electrode;  Q  is  the  charge  per  unit  volume  of  the 
electrode;  Qnvdx  is  the  maximum  charge  that  can  be  ex¬ 
tracted  from  the  electrode. 

The  overall  reaction-rate  at  the  negative  electrode  is 
given  by  Eq.  (1),  which  takes  account  of  the  solid-state 
reactions.  The  overpotential  at  this  electrode  is  defined  as: 

v=4>\  -  <f>2-  (5) 

The  governing  equations  to  develop  the  mathematical 
model  are  shown  in  Table  1  [8,9]. 
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where  the  limiting  current  density  jlim  affects  significantly 
the  polarization  curves  and  is  determined  from  the  dissolu¬ 
tion  rate  of  lead  sulfate,  the  diffusion  rate  of  lead  ions,  and 
the  precipitation  rate  of  lead  crystals. 

The  total  current  density  i  is  the  sum  of  the  current 
density  in  the  solid  phase  /,  and  the  current  density  in  the 
conducting  liquid  phase  i2,  i.e., 

(2) 


2.2.  Initial  and  boundary  conditions 


The  initial  values  for  electrolyte  concentration  and 
porosity  are  as  follows. 


£  ^ref 
£  ~  ^Pb02,ini 
£  ~  £  Pb ,  ini 


(6) 

(7) 

(8) 


where  Eqs.  (7)  and  (8)  relate  to  the  positive  and  negative 
electrode,  respectively.  The  initial  potential  distribution 
can  be  calculated  from  the  equation  for  the  electrode 
kinetic  reaction. 

The  symmetry  conditions,  at  the  centres  of  the  positive 
and  negative  electrode,  are  used  to  define  the  electrolyte 
concentration,  porosity  change,  and  the  voltage  of  the  solid 
phase.  Because  there  is  no  electrolyte  at  these  boundaries, 
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Table  1 

Governing  equations  for  each  region 


Porosity  variation 
Ohm’s  law  in  solid 
Ohm’s  law  in  liquid 
Mass  balance  of  electrolyte 
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Table  2 

Boundary  conditions  for  each  region 
Centers  of  positive  and  negative  electrode 
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Interface  between  separator  and  negative  electrode 
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Table  3 


Coefficients  and  effective  properties  used  in  model  equations 


Positive  electrode 

Reservoir 

Separator 

Negative  electrode 
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the  current  density  of  the  liquid  phase  is  zero.  At  the  centre 
of  the  positive  electrode,  the  potential  on  the  surface  of  the 
current  collector  is  taken  to  be  zero.  The  electrode  kinetic 


Table  4 

Parameters  used  in  calculations  [8] 


Parameter 

Value 

Positive  electrode 

Half  thickness  of  plate 

0.0875  cm 

Maximum  charge  state 

2620  C  cm-3 

Reaction  parameter  (amaxliol  ref) 

0.073  A  cm-3 

“al 

1.15 

“cl 

0.85 

7i 

0.01 

fi 

1.0 

Lead  dioxide  conductivity 

500  S  cm“ 1 

exl 

1.5 

exml 

0.5 

Negative  electrode 

Half  thickness  of  plate 

0.07  cm 

Maximum  charge  state 

3120  Ccm“3 

Reaction  parameter  (tfmaX4**04,ref) 

0.11  A  cm“3 

Q?a4 

1.55 

“c4 

0.45 

74 

0.01 

a 

1.0 

Lead  conductivity 

4.8  X104  S  cm-1 

ex4 

1.5 

exm4 

0.5 

Jlim 

—  102  A  cm-3 

Reservoir 

Thickness  of  reservoir 

0.07  cm 

Separator 

Thickness  of  separator 

0.022  cm 

Porosity 

0.60 

ex3 

1.50 

Electrolyte 

Acid  concentration 

4.9  X  10  3  mol  cm  3 

Transference  number 

0.72 

Partial  molar  volume  of  acid 

45  cm~3  mol 

Diffusion  coefficient 

D  =  exp[(2174.0/298.15) 
-(2174.0/  D]X(1.75 
+  260.0c)X10“5 

Conductivity 

k  =  c  X  exp{l.l  104 

+  (199.475- 16097.781c)c 
+  [(39 16.95 -99406.0c 
-(721860/  T))/T]) 

equation  is  used  to  calculate  the  potential  at  the  negative 
electrode. 

Both  the  flux  of  the  electrolyte  and  the  current  density 
in  the  liquid  phase  are  continuous  at  the  positive  electrode 
solidus  reservoir,  reservoir  solidus  separator,  and  separator 
solidus  negative  electrode  interfaces.  From  the  assumption 
of  electroneutrality,  the  charge  is  conserved  at  these  re¬ 
gions.  The  variation  of  the  porosity,  at  the  positive  elec¬ 
trode  solidus  reservoir  and  separator  solidus  negative  elec¬ 
trode  interfaces,  is  proportional  to  the  rate  of  charge 
reaction.  All  the  current  flows  through  the  liquid  phase, 
because  there  is  no  solid  electrode  at  the  interface  between 
the  reservoir  and  the  separator.  The  boundary  conditions  to 
solve  the  governing  equations  are  shown  in  Table  2  and 
the  coefficients  used  in  the  development  of  the  model  are 
presented  in  Table  3. 

2.3.  Numerical  procedure 

The  space  derivatives  are  discretized  by  the  method  of 
finite  differences  and  the  time  derivatives  are  formulated 
by  means  of  the  Crank-Nicolson  method.  The  non-linear, 
multi-region  problems  are  solved  by  the  Newton-Raphson 
iterative  method  [12]  and  MB  AND  [13].  The  material 
properties  and  the  cell  parameters  used  in  the  simulation 
are  presented  in  Table  4. 


3.  Results  and  discussion 

The  effect  of  various  parameters,  such  as  the  charge-re- 
action  concentration  exponent,  morphology  parameter  and 
limiting  current  density  on  the  charge  performance  is 
investigated  by  using  the  mathematical  model.  The  param¬ 
eter  y  denotes  the  concentration  dependence  on  the  charge 
reaction,  and  £  is  the  morphology  parameter  which  is  used 
to  account  for  the  way  lead  sulfate  covers  the  electrode 
surface.  The  morphology  parameter,  £,  has  a  large  value  if 
the  lead  sulfate  is  spread  well  over  the  electrode  surface. 
The  limiting  current  density,  jlim,  denotes  the  solid-state 
reaction  rate  when  the  lead  precipitates  from  lead  sulfate  in 
the  negative  electrode.  For  constant-current  charge,  the 
charging  current  density  was  changed  and  the  cell  perfor¬ 
mance  was  evaluated.  The  battery  used  in  the  simulation 
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was  discharged  with  1.7  mA  cm-2  for  8  h  at  25°C 
followed  by  a  rest  period  for  1  h  at  the  same  temperature. 

The  effect  of  the  concentration  exponent  of  charge 
reaction,  y,  on  the  cell  voltage  when  the  battery  was 
charged  at  10  mA  cm  is  shown  in  Fig.  1.  The  parame¬ 
ter,  y,  which  corresponds  to  the  reaction-rate  order  in  the 
chemical  reaction,  was  incorporated  to  accommodate  the 
electrolyte  concentration  dependence  on  the  exchange  cur¬ 
rent  density.  The  effect  of  y1  in  the  positive  electrode  is 
simulated  in  Fig.  1(a).  As  yj  is  increased,  the  concentra¬ 
tion  dependence  on  the  electrode  reaction-rate  is  increased 
and  the  voltage  penalty  becomes  higher  at  the  beginning  of 


charge.  After  the  charge  process  has  progressed  at  some 
degree,  it  can  be  seen  that  the  cell  voltage  profiles  at  the 
three  values  of  yx  become  similar.  At  the  initial  state-of- 
charge,  the  electrolyte  concentration  does  not  cause  signifi¬ 
cant  change  of  the  reaction  rate  in  the  positive  electrode. 
By  contrast,  the  electrolyte  concentration  is  an  important 
parameter  in  the  negative  electrode,  where  lead  sulfate 
reacts  with  hydrogen  ions.  The  voltage  penalty  increases 
with  increasing  the  dependence  of  the  electrolyte  concen¬ 
tration  at  the  start-of-charge,  as  shown  in  Fig.  1(b).  The 
slope  of  voltage  increment  becomes  higher  as  y4  is  de¬ 
creased  in  the  negative  electrode. 


(a) 


(b) 


Fig.  1.  Effect  of  charge-reaction  concentration  exponent  for  10  mA  cm  2  charge  at  25°C:  (a)  positive  electrode,  (b)  negative  electrode. 
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The  effect  of  the  morphology  parameter  on  the  cell 
voltage  is  shown  in  Fig.  2.  The  slope  of  the  voltage 
increment  at  the  initial  state  increases  with  increasing 
value  of  in  the  positive  electrode.  This  means  that  the 
electro-active  area  at  the  beginning  of  charge  is  low  and, 
thereby,  increases  the  internal  resistance  of  the  cell.  There¬ 
fore,  the  morphology  parameter  of  the  positive  electrode  t)X 
has  a  strong  influence  on  the  increment  of  the  cell  voltage 
at  the  beginning  of  charge. 

The  effect  of  the  morphology  parameter  in  the  negative 
electrode  on  the  cell  voltage  is  analogous  to  that  in  the 
positive  electrode.  The  higher  value  of  £4  will  have  a 
greater  effect  on  the  voltage  increment  at  the  initial  state- 


of-charge.  This  is  because  nonconductor  lead  sulfate 
spreads  uniformly  over  the  surface  of  the  negative  elec¬ 
trode  and  the  negative  electrode  of  higher  value  of  £4  has 
lower  electrical  conductivity.  When  £4  has  a  low  value, 
the  lead  sulfate  crystal  grows  in  a  needle-like  shape  on  the 
surface  of  the  lead  electrode  during  discharge  [14].  Conse¬ 
quently,  there  is  no  discernible  difference  in  the  electrical 
conductivity  of  lead  electrode  between  the  initial  and  the 
end  of  discharge.  The  increment  of  voltage  is,  therefore, 
much  lower  in  the  case  of  a  high  value  of  £4  than  in  that 
of  a  low  value. 

The  effect  of  the  solid-state  reaction  in  the  lead  elec¬ 
trode  (i.e.,  the  dissolution  of  lead  sulfate,  the  diffusion  of 


(a) 


(b) 


Time  (sec) 


Fig.  2.  Effect  of  morphology  parameter  for  10  mA  cm 


charge  at  25°C:  (a)  positive  electrode,  (b)  negative  electrode. 
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lead  to  the  active  surface,  and  the  precipitation  of  lead 
crystal)  on  the  charge  performance  is  presented  in  Fig.  3. 
The  limiting  current  density  is  determined  by  the  limiting 
dissolution  rate  of  lead  sulfate  or  by  the  limiting  diffusion 
rate  of  lead  ions,  or  by  a  combination  of  both.  For  a  lower 
absolute  value  of  the  limiting  current  density,  the  solid-re- 
action  rate  controls  the  overall  reaction  rate  in  the  negative 
electrode  during  charge  [6].  If  the  solid-state  reaction  rate 
is  much  higher  than  the  charge-transfer  reaction  rate,  the 
absolute  value  of  jhm  becomes  higher.  This  implies  that 
the  predicted  voltage  profile  for  the  solid-state  reaction  at 
the  negative  electrode  will  have  the  same  behaviour  as  in 
the  case  of  including  only  the  charge-transfer  reaction.  The 
cell  voltage  increases  rapidly  with  decreasing  —  jlim  due  to 
the  higher  initial  voltage  during  charge.  Therefore,  if  the 
solid-state  reaction  rate  in  the  negative  electrode  is  low, 
the  charge  efficiency  becomes  lower  due  to  the  high 
internal  resistance  of  the  cell. 

The  distribution  of  the  electrolyte  concentration  in  the 
cell  for  various  limiting  current  densities,  under  the  same 
conditions  as  in  Fig.  3,  is  given  in  Fig.  4.  The  predicted 
profile,  the  electrolyte  concentration  or  a  limiting  current 
density,  jhm,  of  —  1.0  A  cm-3  is  shown  in  Fig.  4(a).  If  the 
absolute  value  of  the  limiting  current  density  is  low,  the 
rate  of  solid-state  reaction  controls  that  of  the  electrode 
reaction  in  the  negative  electrode.  It  is  found  that  only 


16%  of  the  previous  discharge  has  been  restored.  Although 
the  concentration  gradient  in  the  negative  electrode  is  not 
steeper,  the  polarization  resistance  becomes  higher  due  to 
the  slow  rate  of  the  solid-state  reaction.  In  the  electric 
field,  the  electrolyte  concentration  in  the  positive  electrode 
is  higher  than  that  in  the  negative  due  to  the  migration 
resistance  of  the  HSO^  ions  [15].  The  limiting  current 
density  is  taken  to  be  —  5.0  A  cm-3  in  Fig.  4(b).  Thus,  the 
solid-state  reaction  rate  is  faster  in  Fig.  4(b)  than  in  4(a). 
At  100  s,  for  the  electrolyte  distribution  there  is  no  signifi¬ 
cant  difference  between  cases  (a)  and  (b).  The  charge 
returned  is  about  74%.  When  the  limiting  current  density  is 
set  to  —  10.0  A  cm-3,  as  shown  in  Fig.  4(c),  the  behaviour 
of  the  electrolyte  concentration  is  similar  to  Fig.  4(b)  and 
about  79%  has  been  returned.  If  the  solid-state  reaction 
rate  is  very  fast,  the  rate-determining  step  in  the  electrode 
reaction  becomes  the  charge-transfer  rate,  as  represented  in 
Fig.  4(d).  It  is  shown,  in  this  case,  that  the  charge  restored 
is  about  84%.  From  the  results  given  in  Fig.  4(a)-(d),  it  is 
concluded  that  the  limiting  current  density  affects  signifi¬ 
cantly  the  initial  voltage  rise,  while  it  has  only  a  very  little 
influence  on  the  concentration  gradient  in  the  cell,  when 
the  charging  current  density  is  constant. 

The  cell  voltage  is  demonstrated  in  Fig.  5  for  charging 
current  densities.  When  the  charging  current  densities  are 
set  to  1,  5,  10  and  20  mA  cm  ,  respectively,  the  times  to 


Concentration  (mol/l)  B  Concentration  (mol/l) 
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Fig.  4  ( continued ). 
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voltage  cut-off  are  47,762,  8295,  3255  and  424  s,  respec¬ 
tively.  As  the  charging  current  is  increased,  the  internal 
resistance  of  the  cell  is  higher  due  to  the  limitation  of  mass 
transfer  of  lead  ions  [16].  With  time,  therefore,  high-rate 
charging  shows  a  steeper  rise  in  the  cell  voltage.  This  is 
because  the  internal  resistance  has  a  greater  effect  on  the 
voltage  penalty. 

The  concentration  distribution  in  the  cell  is  demon¬ 
strated  with  the  same  conditions  as  in  Fig.  5.  When  the 
current  density  is  1  mA  cm-2,  the  distribution  of  elec¬ 
trolyte  concentration  with  time  is  represented  in  Fig.  6(a). 
The  charge  returned  is  about  98%  of  the  previous  dis¬ 
charge.  The  concentration  in  the  cell  during  charge  was 
relatively  uniform,  which  indicates  the  lower  internal  resis¬ 
tance  of  the  cell.  The  limiting  current  density  is  —3.0  A 
cm-3  and  the  solid-state  reaction,  in  this  case,  is  not 
ignored  as  shown  in  Fig.  3.  After  the  end  of  the  charge 
process,  the  concentration  difference  between  the  end  of 
the  positive  electrode  and  that  of  the  negative  electrode  is 
very  low;  thus,  the  increment  of  the  internal  resistance  due 
to  the  concentration  gradient  is  not  larger.  When  the 
charging  current  density  is  5  mA  cm-2,  about  85%  of  the 


discharge  is  restored,  as  shown  in  Fig.  6(b).  As  the  charg¬ 
ing  current  density  increases,  the  time  to  the  voltage 
cut-off  is  shortened  but  the  concentration  gradient  is  in¬ 
creased.  In  particular,  the  concentration  gradient  is  maxi¬ 
mum  at  the  interface  between  the  positive  electrode  and 
the  reservoir.  For  a  charging  current  density  of  10  mA 
cm  ,  approximately  67%  of  the  previous  discharge  is 
returned.  The  concentration  distribution  in  the  positive 
electrode  shows  that  the  gradient  is  much  steeper  than  that 
in  the  other  regions.  As  shown  in  Fig.  6(b),  the  concentra¬ 
tion  gradient  is  much  higher  at  the  interface  between  the 
positive  electrode  and  the  reservoir.  Therefore,  the  resis¬ 
tance  of  the  mass  transfer  in  this  region  increases  with 
increasing  charging  current  density.  On  charging  the  bat¬ 
tery  at  high  rates,  the  charging  efficiency  is  lower  due  to 
the  concentration  polarization  and  the  voltage  rise  at  the 
beginning  of  charge.  At  a  charging  current  density  of  20 
mA  cm-2,  the  concentration  profile  in  the  cell  has  a 
nonuniform  distribution,  as  demonstrated  in  Fig.  6(d).  Note 
that  the  concentration  near  the  current-collector  of  the 
positive  electrode  is  markedly  higher  and  there  is  a  site  at 
which  the  electrolyte  concentration  is  largely  unchanged  in 


Time  (sec) 

Fig.  5.  Effect  of  charging  current  density  for  ylim  = 


-3.0  A  cm-3  at  25°C. 
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Fig.  6.  Predicted  profiles  of  acid  concentration  for  j'lim  =  —  3.0  A  cm-3  at  25°C:  (a)  1  mA  cm-2,  (b)  5  mA  cm-2,  (c)  10  mA  cm-2,  (d)  20  mA  cm 


26 


S.C.  Kim,  W.H.  Hong  /  Journal  of  Power  Sources  89  (2000)  15-28 


Distance  (cm) 

Fig.  6  ( continued ). 


S.C.  Kim,  W.H.  Hong  /  Journal  of  Power  Sources  89  (2000)  15-28 


27 


the  reservoir.  This  is  because  the  transfer  resistance  of 
charge  produced  both  in  the  positive  and  negative  elec¬ 
trode  is  higher. 


4.  Conclusions 

The  effect  of  various  parameters  on  the  cell  perfor¬ 
mance  for  constant-current  charge  is  investigated  by  means 
of  a  mathematical  model.  The  results  from  applying  the 
model  can  be  summarized  as  follows. 

(i)  A  higher  value  of  the  concentration  exponent  of 
charge  reaction,  y,  in  the  negative  electrode  has  a  greater 
effect  on  the  charge  performance  than  that  in  the  positive 
electrode,  due  to  increase  in  the  voltage  penalty. 

(ii)  The  slope  of  the  voltage  rise  at  the  initial  state 
increases  with  increasing  value  of  £,  the  morphology 
parameter.  This  means  that  the  active  surface  area  at  the 
beginning  of  charge  is  low,  and  causes  an  increase  in  the 
internal  resistance  of  the  cell. 

(iii)  The  charge  at  a  low  value  of  —  jlim  is  inefficient 
because  the  charge  performance  is  limited  by  the  solid-state 
reaction  rate  in  the  negative  electrode. 

(iv)  As  the  charging  current  density  is  increased,  the 
concentration  gradient  has  a  maximum  value  at  the  inter¬ 
face  between  the  positive  electrode  and  the  reservoir. 
Thus,  the  internal  resistance  of  the  cell  becomes  higher  due 
to  the  limitation  of  the  transfer  rate  of  the  electrolyte. 


5.  Symbols 

a  active  surface  area  of  electrode  (cm2  cm-3) 

ax,  a2  coefficients  used  in  governing  equations  and 
boundary  conditions 

amax  maximum  specific  active  surface  area  of  elec¬ 

trode  (cm2  cm-3) 

c  concentration  of  binary  electrolyte  (mol  cm-3) 

cref  reference  concentration  of  the  binary  electrolyte 

(mol  cm-3) 

D  diffusion  coefficient  of  the  binary  electrolyte 

(cm2  s-1) 

F  Faraday’s  constant  (96,487  C  mol-1) 

i  total  applied  current  density  based  on  projected 

electrode  area  (A  cm-2) 

i}  current  density  in  solid  phase  (A  cm-2) 

i2  current  density  in  conducting  liquid  phase  (A 

cm-2) 

zol  ref  exchange  current  density  at  cref  for  positive 

electrode  (A  cm-2) 

/o4ref  exchange  current  density  at  cref  for  negative 

electrode  (A  cm-2) 

j  reaction  current  per  unit  volume  of  electrode  (A 

cm-3) 

yjim  limiting  current  density  for  negative  electrode  (A 
cm-3) 


MWi  molecular  weight  of  species  i  (g  mol-1) 

Q  charge  density  in  electrode  (C  cm-3) 

<2max  theoretical  maximum  capacity  (C  cm-3) 

R  universal  gas  constant  (8.3143  J  mol-1  K-1) 

t  time  (s) 

transference  number  of  H+  with  respect  to  sol¬ 
vent  velocity 

T  absolute  temperature  (K) 

f/pbo2  equilibrium  potential  at  cref  for  positive  elec¬ 

trode  (V) 

v  *  volume  —  average  velocity  (cm  s  - 1 ) 

v  distance  from  centre  of  positive  electrode  (cm) 

Greek  letters 

aal  anodic  transfer  coefficient  for  positive  electrode 
acl  cathodic  transfer  coefficient  for  positive  elec¬ 

trode 

aa4  anodic  transfer  coefficient  for  negative  electrode 

ac 4  cathodic  transfer  coefficient  for  negative  elec¬ 

trode 

y1  concentration  exponent  for  positive  electrode 

y4  concentration  exponent  for  negative  electrode 

s  porosity 

£sep  porosity  of  separator 

^Pbo,  porosity  of  positive  electrode  at  initial  state-of- 

charge  reaction 

^Pb,ini  porosity  of  negative  electrode  at  initial  state-of- 
charge  reaction 

C1  morphology  parameter  for  positive  electrode 

C4  morphology  parameter  for  negative  electrode 

r]  total  local  overpotential  with  respect  to  equilib¬ 

rium  potential 

k  electrolyte  conductivity  (S  cm-1) 

Pi  density  of  species  i  (g  cm-3) 

(T j  conductivity  of  electrode  matrix  (S  cm-1) 

<pl  potential  in  electrode  matrix  (V) 

<j)2  potential  in  solution  (V) 

Superscripts 

ex  exponent  on  porosity 

exm  empirically  determined  constant  for  tortuosity  of 

solid  matrix 
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